Non-uniform glassy electronic phases from competing local orders 
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We study non-uniform states and possible glassiness triggered by a competition between distinct 
local orders in disorder free systems. Both in Ginzburg-Landau theories and in simple field the- 
ories, such inhomogeneous states arise from negative gradient terms between the competing order 
parameters. We discuss applications of these ideas to a variety of strongly correlated systems. 
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I. INTRODUCTION. 

Accumulated experimental evidence strongly suggests 
that in many correlated electronic systems, different 
types of ordering phenomena compete and coexist over 
a wide range of tunable parameters. The most ubiqui- 
tous such cohabitation is between magnetic and super- 
conducting orders. Itinerant antiferromagnctism (AFM) 
coexists with superconductivity in the 115 heavy fermion 
series (CeMIus, where M =Co,Ir, or In)i. In UPt^, su- 
perconductivity emerges at Tc ~ 0.5K from a strongly 
correlated heavy electron state with small moment AFM 
below 6K— . In some of the high-T^ cuprates charge den- 
sity order coexists with spin density order ( "stripes"^!^) 
and may be relevant to the onset of the superconductivity 
and quantum critical behavios^. Recent measurements 
indicate that in URu2Si2 there is a proliferation of com- 
peting phases under an applied magnetic field^. Exper- 
iments also suggest multiple phases in the skutterudite 
superconductor PrOs4Sbi2^, manganitesi^, and a num- 
ber of other materials. Two trends are common to these 
experimental findings. First, the coexistence of different 
orders is often inhomogeneous. Second, this coexistence 
is frequently most pronounced near a Quantum Critical 
Point (QCP), where the transition temperature for one 
of the order parameters vanishe o^d^ . 

Additionally, dynamics of compounds with inho- 
mogeneous coexistence of distinct orders is often 
glass^ r^^d'^d^d^ . In some systems, such as manganitesi^ 
the glassy behavior is, most likely, due to disorder upon 
doping the system. In others, including cuprates, the 
glassiness may be self-generated (not simply due to dop- 
ing disorder—), and arise out of competing interactions 
at different length scalesi^. The question remains, how- 
ever, whether inhomogeneous and/or glassy behavior can 
arise out of a theory with local interactions and no disor- 
der. In this article we address this question for a class of 
Ginzburg-Landau theories with competing order param- 
eters. A comprehensive survey of classical systems with 
frustration and no disorder that display glassy behavior 
and proliferation of inhomogeneous ground states, can be 
found in Refs.IHi^. 

We study a minimal Ginzburg-Landau (GL)^ theory 
which includes amplitude-gradient coupling between two 



distinct local orders, and find the conditions for resul- 
tant inhomogeneous phases. A related interesting work 
examining gradient coupling in GL theoriesS^ appeared 
slightly after the initial dissemination of our results^. 
Extensions of the GL gradient couplings considered here 
are found in some studies of the supersolid transitional. 
We show that, for a range of parameters described be- 
low, our theory maps onto an effective model that is likely 
to exhibit glassiness. Whether a particular system does 
or docs not show glassy behavior upon cooling depends 
on the rate of temperature change and other dynamical 
variables that are not part of our equilibrium analysis. 
However, our approach allows us to conclude whether a 
glassy phase is possible and likely to occur. In this we 
follow the established approaches in the fieldi^. 

The mapping that strongly suggests glassiness in our 
approach is to a Brazovskii-like model for one of the order 
parameters. The Brazovskii model^ for a single compo- 
nent order parameter is defined by a GL functional of the 
form 
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in momentum (fc) space with V the volume of the system. 
In Eq.([T|), the ellipses denote cubic, quartic, and higher 
order terms in the order parameter field $. As the mass 
term, rg changes sign, the transition to a broken sym- 
metry state $ ^ 0, involves the appearance of structures 
characterized by a finite wavenumbcr on a shell of radius 
g > 0). Structures that satisfy definite commcnsura- 
bility relations amongst the wavcnumbers arc most pre- 
ferred. In Ref. [25I Brazovskii found that large phase space 
available for fluctuations around the minimizing shell al- 
ters the character of the transition to the ordered state 
once the fluctuations are accounted for, and suggested 
that it becomes first order. Thermal fluctuations renor- 
malizc the cubic terms of the GL theory. More recent 
replica calculation a^^'^^i^^'^^ showed that the model has 
extensive configurational entropy, indicating proliferation 
of modulated low-energy states, and strongly suggesting 
slow dynamics and glassiness under generic experimental 
conditions. Once again, these replica calculations only 
establish that glassiness is a plausible and likely alterna- 
tive to the first order transition into a uniformly mod- 
ulated phase. Whether a finite temperature Brazovskii 
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transition does or does not transpire before the system 
undergoes a dynamical arrest (the glass transition out- 
lined below) depends on microscopic details of the model. 
The known theoretical techniques (SCSA, DMFT, and 
others) do not enable the proof of a glassy phase. These 
methods only enable us to determine whether a glassy 
phase is possibl e 

Below we find the mapping of systems with compet- 
ing orders to Brazovskii type models. This mapping 
allows us to (i) Find resultant inhomogeneous phases 
in the GL analysis; (ii) Include fluctuations via a self- 
consistent field theory to establish that one of the two 
scenarios is realized: (a) the critical temperature for the 
onset of non-uniform states is suppressed to zero, sug- 
gesting that these states are more likely to be observed 
near a QCP; or, alternatively, (b) fluctuations lead to a 
low temperature Brazovskii transition; (iii) Appeal to ex- 
isting replica calculation results to confirm the extensive 
configurational entropy associated with these incommen- 
surate structures in disorder free systems with competing 
local orders, which strongly suggests slow dynamics and 
glassiness. Finally, wc comment on possible realizations 
of our model and applicability of the results to itinerant 
electronic systems. 



II. GINZBURG-LANDAU THEORY: 
INSTABILITY OF UNIFORM COEXISTENCE. 

To empirically account for competing orders, we ana- 
lyze the Ginzburg-Landau (GL) functional with two or- 
der parameters, <I>i and $2, which we will choose to be 
real and scalar without loss of generality. We remark that 
our very general GL approach applies to various types of 
order parameters. Of course, the symmetry, number of 
components of the order parameters etc. changes. Nev- 
ertheless, the conclusions arc generally much the same. 
The uniform part of the free energy is J-q = J c?x_fo where 



^^o = ^i$ir + ^i<i>2i^ + |i$ii 



l$2 



-fH>2l 



In the spirit of the GL theory, ri_2 — 0-1.2(2^ ~ Ti, 2), with 
Ti the mean field transition temperatures. All other co- 
efficients are taken to be temperature- independent. The 
quadratic coupling of the order parameters is allowed for 
all symmetries. We consider competing orders, i > 0, so 
that the uniform coexistence region (<l>i 7^ 0, $2 7^ 0) oc- 
curs only below the lower of the transition temperatures, 
and for u > t^. In that case, the values of the fields mini- 
mizing the free energy are = (r2t — riu)/ (u — t^), and 
^"2 ~ {fit ~ ''"2) I {u ~ t^). In disorder free systems, the 
only alternative to the uniform coexistence is phase sep- 
aration unless non-trivial gradient terms arc preseni*2S. 
Therefore we include the inhomogeneous contribution to 
the free energy, J^g = J (ixFg(x), where 



.(3) 



Here, we included general symmetry allowed low order 
gradient terms. To flesh out the quintessential physics in 
what follows, we set 511 = 1722 = .921 = 0,1712 > 0, and 
pi = 0. This is the essential aspect of the model that 
allows us to investigate the appearance of the inhomoge- 
neous states. The coupling of the form — .gi2|'I'iP|V<f>2p 
implies that in the effective theory for the order param- 
eter $2 the coefficient of the gradient term, 1 — gi2|'I'iP, 
may become negative, making the transition of the Bra- 
zovskii type. We now investigate when this is possible. 

With F(x) = Fo + Fq, the order parameter profiles sat- 
isfy the Euler-Lagrange equations, [V • (9F/9(V$i))] = 
(dF/d^i). By constructing inhomogeneous variational 
states whose free energy is lower than the minimum 
amongst all possible uniform configurations, we prove 
that the uniform solution is unstable towards the appear- 
ance of inhomogeneities. Wc study the phase diagram of 
the model assuming that the mean field transition tem- 
peratures Ti can be tuned by an external parameter, x 
(pressure, doping, magnetic field etc.), as shown in Figll] 
with Ti{x) monotonically decreasing, and T2{x) mono- 
tonically increasing. That is. 



Ti = T^^ - aix, 



(0) 



a2X, 



(4) 



with ^i^d ai^2 positive constants. 

We first concentrate on the region Ti > T2. Upon 
lowering the temperature, the first transition is into the 
uniform state with $2 = and ^•i(x) ~ — ri. Conse- 
quently, below the Tq ~ Ti — l/(gi2ai) the coefficient 
of the |V$2p term becomes negative indicating the ten- 
dency towards the development of an inhomogeneous $2 
phase. The structure of this modulation depends on the 
difference Tq — T2. If this difference is sufficiently large, 
it is disadvantageous to create non-vanishing bulk aver- 
age of $2 • Local "bubbles" of the order may appear upon 
lowering T , but their study is not our focus in the present 
work. 

.(2) In order to make the connection with the slow dy- 
namics and Brazovskii transition, we study the on- 
set of the periodically modulated phase of the form 
<i>2(x) = 92C0s(qi • x). Of the numerous contending 
low (free) energy configurations, wc will focus on ana- 
lytically tractable modulated structures; we do so in or- 
der to obtain stringent variational bounds that we are 
able to extremize, and based on the original analysis 
that showed the single modulation structures are most 
advantageous^^. In the regime Tq > T2 minimization of 
the GL functional with respect to both q and 02 gives 
the transition temperature 



Tc2 = Tq ~ (.gi2ai) 
_ pa2 - tpai 



z^+'^ + 2pa2{Tq-T2)-z 

912 



9120.1 



(5) 



to the phase O2 ^ with modulations at a finite wave 
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FIG. 1: The phase diagram obtained from the Ginzburg- 
Landau expansion. The lines T\^2 denote the bare mean field 
transition temperatures as a function of a tuning parameter 
X. Tq is defined in text. An inhomogeneous phase appears 
below Tc2- Double line denotes the first order transition. 



vector, 
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In the regime Tq < Ti^2 the first transition, at Tg < T, 
occurs into a spatially homogeneous phase. We next in- 
vestigate the phase diagram for the more general varia- 
tional ansatz <i>2°'"(x) = $2 + 62 cos(qi ■ x). Introduction 
of spatial modulations reduces the condensation energy 
and therefore is unfavorable, unless compensated by a 
significant gain due to the negative gradient term. As 
a result we find a (generically first order, but dependent 
on the magnitude of the coefficients in the GL expansion) 
transition from the homogeneous to modulated, with a 
finite q, phase at low T. In Fig.lJT]), we show the phase 
diagram of Eqs.®-® for t ~ 0. Of course, since we 
allowed only for the restricted variational states in the 
above analysis, our bounds are more potent for the global 
free energy minima - $2 is strictly inhomogeneous for 
all T < Tc2{x); unrestricted inhomogeneous states (not 
bound to the form of <I>™'") may extend to temperatures 
somewhat higher than Tc2{x). 



III. SELF-CONSISTENT FIELD THEORY FOR 
COMPETING ORDER PARAMETERS. 



To improve on the GL analysis and incorporate the 
effect of fluctuations self-consistently, we generalize our 
model to 71-component vector fields and utilize a large n 
expansion. As well known, the n = 00 limit is equiva- 
lent to the spherical model describing single component 
(scalar) particles^. The physical engine for the inhomo- 
geneities is, as in preceding section, the amplitude gra- 
dient coupling which drives non-uniformities in $2 once 
$1 is finite. For a finite ^'i(x) = $1, the effective free 



energy for 4>2 is 

d'^k \j2 



e//;2 



(v + l5*?) + (l-5i2$?)fc'+pfc1 



X $2(k)$2(-k) 
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d'^ki d'^k2 d'^ks 



4 J {2t:Y {2ttY {2t:Y 
X 4>2(ki)$2(k2)$2(k3)$2(-kl -k2 -kg), (7) 

where d is the dimensionality of the system. The bare 
inverse Green's functions are given by G^"^ = [^2/2 -|- 
t/2-\-{l—gi2^\)k^+pk'^]. Incorporating fiuctuations self- 
consistently, we have = [r2/2 + (1 - gi2^\)k'^ +pk'^] 
where, by the Dyson equation, r2/2 = r2/2 -I- t/2 + S. 
To lowest order in 1/n, the self-energy is given by YP = 
J ■^I^G(k), see Ref. [S^ This leads to a self-consistency 
equation for 7^2- Similar self-consistency equations appear 
for $1; before the transition to an ordered $2 state, <I>^ = 
—ri. A phase transition to an ordered state $2 
occurs when the Green's function acquires a pole on the 
real k axis. If the pole is at fcmin = 0, the transition is 
to a uniform phase of <i>2; if the pole first appears for 
^min 7^ 0; the transition is into a modulated phase. 

When [1 — 512 $1] > the minimum of is always 
at fc = 0, and both $1 and <I>2 may exhibit uniform or- 
ders. On the other hand, if [1 — gi2't>i] < 0, the mini- 
mum for the $2 inverse Green's function, G'^^(k) occurs 
at kjnin = ^[1 ~ 5i2'5'i]/(2p) leading to a real axis pole 
when r2 = r2 min = [1 - gi2^'f]^/(2p2)- The quartic G~^ 
has two pairs of complex conjugate poles in the k plane 
which lie on a circle of radius p = (r2/(2pi2))^/'*. The 
finite real component of the poles means that the corre- 
lation function ($2(x)<I'2(y)} exhibits sinusoidal modula- 
tions in addition to exponential decay. The modulation 
and correlation lengths are given, respectively, by 



h = 47r[V?^ + (1 - .gi2$?)/2]-i/2, 

6 = 2[v^-(l-.gl2<i>?)/2]-l/^ 



(8) 



with $1 the uniform competing order field. 

Irrespective of the spatial dimensionality, whenever 
[1 — 512 $1] < 0, as r2 r2 min the self-energy diverges 
as S ^ (r2 — r2 min)""'^^^- The phase transition which 
would occur (at the mean field level) when r2 = f2 min 
is thwarted by the divergence of the self energy due to 
fluctuations. This implies that Tc = similar to systems 
with competing long range interactions^. However, fi- 
nite n corrections (especially for n = \), may make the 
transition temperature finite. In this case, this low tem- 
perature transition for <i>2 is of Brazovskii type, with a 
shell of minimizing modes. 

A similar analysis holds for competing local orders in 
large n quantum systems by extendin g^^'^'^ . For bosonic 
fields, after a summation over Matsubara frequencies, the 
$2 correlator is 



G2(k) = 



[f + (l-.9l2$?)fc2+pfc4]/(fcBr) 



f + (l-gi2$?)fc2+pfc4 
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with nsix) = [cxp(2:) — 1]""'^. Order, at large n, is still in- 
hibited in the quantum rendition of our system, although 
the divergence of the self-energy in this case is less severe 
than for its classical counterpart. In the bosonic system, 
due to integration over imaginary time, the <I>2 self energy 
diverges as —\n\r2~r2 mini when ¥2 r2 min), whereas 



in the classical system it diverges as [r2 — r2 



-1/2^ 

Thermal fluctuations in classical n = 2 systems (e.g. 
complex scalar fields) also lead to a divergence of the 
-ln|r2 - r2 mini] typei^i^^i^ The order probed by $2 
is stabilized if the degeneracy of the minimizing wave- 
numbers is lifted by augmenting the rotationally symmet- 
ric Hamiltonian by additional lattice point group symme- 
try terms. In such instances, the critical temperature of 
the large n system often remains anomalously low, and 
attains its minimal value exactly at the onset of incom- 
mensurate order (e.g., when the minimizing modes are 
of vanishing norm((7 — > 0+))^. Degeneracy can also be 
lifted by an external field which lowers the full rotational 
symmetry of H to a lower rotational symmetry in a plane 
orthogonal to the field direction. For small n, a "quan- 
tum" finite temperature Brazovskii transition may occur 
for $2. 



IV. SLOW DYNAMICS AND CLASSINESS. 

The key conclusion of the previous section was that 
our model of competing orders, Eqs.([2l)-(l3|), maps onto a 
Brazovskii-like model for the subdominant order param- 
eter. The transition temperature for the onset of this 
order is then suppressed relative to the mean field con- 
clusions due to a large phase space available to low energy 
thermal and quantum fluctuations. There are two pos- 
sibilities for the transition itself. It may take place as 
a first order fluctuation-induced transition as originally 
envisioned by Brazovskii. Recent work on a model equiv- 
alent to our Eq. ([7]) suggests that a glass transition may 
be realized as an alternativei^i^. 

When applied to Eq.Q, the self-consistent screening 
approximation shows^i^ that the configurational en- 
tropy, Sc ~ ks^ogNm, with Nni the number of the 
metastable states, is extensive (proportional to the vol- 
ume) over a finite temperature range (Ta > T > Tk), 
which depends on the coefficients of our GL expansion. 
This entropy is due exclusively to the inhomogeneous 
field $2(x) triggered by the competing uniform order 
<I>i(x) = $1. At the onset (T = Ta) 



Sc{Ta) « CfcB(,9l2$l - l)^V, 



(9) 



where, in 3D, the numerical constant C ~ 1.18 x 10""^ 
and V is the volumei^. Our resulting effective model for 
the inhomogeneous order is exactly equivalent to that 
of Ref. |2^, after shift in the term and recognition 
that the self-consistently determined effective tempera- 
ture (r), depends on the dominant order $1. The reason 
for equivalence is that in the presence of only quartic and 
biquadratic terms, the propagator lines for the fields $1 




FIG. 2: Schematic phase diagram beyond the GL theory. 
Here we highlight the possibility of a glassy phase triggered 
by the competition of two local orders. Alternatives include 
the first order Brazovskii transition into a modulated state, 
or a transition with a severely suppressed Tc- 



and $2 are continuous and only allow for self energies 
with the same field index as for a single field problem 
with shifted parameters. In Ref. [28l the dynamical mean 
field theory calculation yielded extensive Sc for the sin- 
gle component problem, and hence precisely the same 
conclusion is applicable to our model. 

The extensive value of Sc implies that N,n oc e^, and 
strongly suggests glassiness for T < Tk>^. The condi- 
tion for possible glassiness formulated in Refs. [l8ll28l is 
that that ratio of the coherence length to the modula- 
tion scale exceeds a number of order two. As seen from 
Eq.® in our model at low temperatures ^2/^2 > 2 sat- 
isfying this condition. Once again, the realization of the 
glassy phase depends on the details of dynamics in a par- 
ticular measurement, but extensive entropy makes such 
an outcome likely. 

The high degree of low temperature entropy can 
be made rigorous. In all large n (and several Ising) 
systemsS^, the extensive configurational entropy found 
at higher temperatures by replica calculations is sup- 
planted by a ground state degeneracy scaling as the sur- 
face area of the system (Sground oc gd-iy(d-i)/d)26^ 
in d spatial dimensions. By explicit construction, these 
systems can be shown to possess a multitude of zero en- 
ergy domain walls^^. These low temperature excitations 
going hand in hand with a multiple metastable low en- 
ergy states. Numerical simulations of single component 
systems in similar clasical models of liquids also report 
exceptionally sluggish dynamics^i^i with strong indica- 
tions of glassiness^. Thus, the non-uniform structures 
arising in our model of competing order parameters natu- 
rally exhibit slow dynamics and is likely to become glassy. 

Summarizing, the field theoretical analysis accounting 
for fluctuations around the inhomogeneous minimizing 
structure extends the GL picture and strongly suggests 
the phase diagram shown in Fig. ([2]). For the low n sys- 
tems of relevance, the low temperature first order Bra- 
zovskii transition can be pre-empted by a transition into 
a glass. 
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V. RELEVANCE TO ELECTRONIC SYSTEMS. 

We showed that, when there is a competition between 
two order parameters of different origin, and when a gen- 
eral symmetry allowed gradient-amplitude coupling in a 
local theory is negative {even if of moderate magnitude), 
the coexistence of two orders is inhomogeneous, and, gen- 
erally, either the dynamics of the system is slow or a 
first order Brazovskii transition occurs. Crucially, even 
though we start with a local theory, the inhomogeneous 
coexistence leads to a low-energy theory of the same class 
as considered in models of self-generated glassiness due 
to competing length scales of interactioni^i^l, although 
the origin of the phenomenon is very different. 

Moreover, the transition temperatures for both order 
parameters are suppressed compared to the mean field 
value. We emphasize that the gradient-amplitude cou- 
pling is required to stabilize an inhomogeneous state: in 
its absence only uniform or phase-separated configura- 
tions are thermodynamically stable, as has been shown 
for stripe orders^S. Therefore in competing coexisting 
phases, Tc is lower than in other parts of the phase di- 
agram, structure factor measurements will indicate non- 
uniform order, and dynamical measurements likely dis- 
play slow dynamics. The natural question to ask is what 
systems offer the best chance for realization of the model 
considered above. 

Cuprates provide one obvious example of such com- 
peting orders when static low temperature spin and 
charge density waves (stripes) are inhibited in the pres- 
ence of superconducting order. In these materials, STM 
measurements^ indicate incommensurate coexistence of 
superconductivity and a pseudo-gap state at nanoscale 
level; however, the dynamics in this situation is strongly 
energy dependent, which suggests that the mapping on a 
simple GL theory with temperature-independent coeffi- 
cients is insufficient. At least in one example the scaling 
form of the dielectric function in the glassy state goes 
smoothly to quantum critical scaling as the glass transi- 
tion temperature tends to zercJ^. 

Heavy fermion systems provide perhaps the best 
chance for observing the phenomena described here. In 
materials of the 115 family proximity or coexistence of 
antiferromagnetic and superconducting phases is now 
well established^!^, and experiments indicate an inho- 
mogeneous coexistence of the two orders in a magnetic 
field^i^. Moreover, there is strong evidence that Cd 
and Hg dopants^^iiS create antiferromagnetic regions in 
their vicinity^i^i^, suggesting that the system is on the 



border of inhomogeneous coexistence of two orders. The 
Neel temperature drops precipitously if superconducting 
transition occurs first^. No dynamical measurements 
have yet been carried out in the relevant regime of the 
phase diagram, but it would be interesting to see if, for 
example, in CeRhlns under pressure the spin dynamics 
as determine by NMR shows signatures of slowing down 
or freezing at low temperatures. 

In several systems the inhomogeneous coexistence was 
proposed in the presence of coupling terms that exist 
only under special circumstances^^i^i^. A particularly 
relevant example are manganites where the coupling due 
to deviation from half-filling that promotes the inhomo- 
geneous coexistence of the magnetic and charge orders 
was proposed recently based on considerations similar to 
ours^. As mentioned above, in these materials glassi- 
ness may emerge due to bona fide disorder, and not be 
self-generated. Non-trivial couplings appear in some of 
the multiferroic materials, e.g. spiral magnets such as 
RMnOs with R= Tb, Ho, Dy.^ 

It is important to note that, if we extend the treat- 
ment to include external parameters such as strain and 
field to act as a massive (i.e. with fluctuations towards 
order but no symmetry breaking since the quadratic coef- 
ficient in the GL expansion remains positive) "competing 
orders" within the GL framework, the resulting inhomo- 
geneous state only occurs for moderately large coupling. 
One candidate for such a scenario is MnSi, where the 
low energy theory exhibiting these features has recently 
been put forward on the basis of Dzyaloshinskii-Moriya 
coupling^. 

In conclusion, we believe that many of the observed low 
temperature transitions, inhomogeneities, and slow dy- 
namics/giassiness found in strongly correlated electronic 
systems are a natural consequence of competing local or- 
ders. As we illustrated, competing local orders may trig- 
ger inhomogeneities with likely first order transition or 
possible glassiness. In our calculations, the proliferation 
of incommensurate ground and metastable states is the 
common origin of both the dramatic lowering of the tran- 
sition temperature or viable first order Brazovskii tran- 
sition and possible glassy dynamics. 
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